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P-i! We study the temperature dependence of quark antiquark correlations in the space hke di- 

P^l rection. In particular, we predict the temperature dependence of space like Bethe-Salpeter 

amplitudes using recent Lattice gauge data for the space like string potential. We also 

investigate the effect of the space like string potential on the screening mass and discuss 

possible corrections which may arise when working with point sources. 
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1 Introduction 

It is generally believed that hadronic matter undergoes a transition to a new phase, the 
so-called quark-gluon-plasma (QGP), which is a gas of weakly interacting quarks and 
gluons.(For a general review see refs. fl], 0). As far as thermodynamic properties are 
concerned this statement seems to be supported by lattice gauge calculations (LGC)[|, 
which show a sharp rise of the energy density at the critical temperature Tg. Above Tc, 
the energy- and entropy density obtained in LGC are compatible with a perturbative gas 
of quarks and gluons. Moreover the quark condensate vanishes above T^ indicating the 
restoration of chiral symmetry. Furthermore, above T^ screening masses extracted from 
spatial point to point correlation functions assume the values of Mg^ = nnT |^, |^, n being 
the number of quarks involved. The (tt, a) correlators give screening masses smaller than 
the free value indicating considerable residual interactions in these channels above Tc. 
This is not too surprising as these modes can be viewed as long range fluctuations in the 
order parameter of a second order chiral restoration transition. 

There exist, however, several phenomena which indicate a nontrivial structure even 
above Tc. Simple thermodynamic considerations 0, for instance, show that about 50% 
of the zero temperature gluon condensate still remains condensed above the transition. 
Moreover, while the large distance behavior of space like point to point correlators is 
consistent with free quarks, the measurement of Bethe-Salpeter amplitudes in the space 
like direction shows very strong correlations between quarks and antiquarks 0. It has 
been suggested |^ and demonstrated explicitly p that these 'observed' correlations can be 
related to the so-called space-like string-tension, which remains finite even above Tc |iy . 



This space like string tension is extracted from the expectation value of a Wilson loop, 
which loops only in the spatial direction. Contrary to the time like Wilson loop, which 
measures the electric interactions and thus can be related to the heavy quark potential, the 
space like Wilson loop is sensitive to magnetic interactions. The associated string tension, 
therefore, does not relate to a heavy quark potential. It rather provides a measure for the 
strength of a current-current type interaction, analogous to that associated with Amperes 
law in electrodynamics. In contradistinction to the time like string tension, which vanishes 
above Tc, the space like string tension remains finite at all temperatures [^, |T^ 

The effect of the spatial string tension on a quark antiquark pair can be best under- 
stood if one transforms to so called 'funny space', where the euclidian time direction and 
one of the spatial directions, say the 'z '-direction, are interchanged. At zero temperature, 
of course, real euclidian space and funny space are identical. At finite temperature, how- 
ever, bosons and fermions follow periodic and antiperiodic boundary conditions in the 
temporal direction, respectively. In funny space these boundary conditions are moved to 
the 2;/-direction with no boundary conditions in the Tf direction. Consequently a sys- 

^For a collection of the most recent LGC results see ref. [|j. 



tern at finite temperature in real space transforms to one at zero temperature in funny 
space where quarks and gluons form standing waves in the ^/-direction with momenta 
Pq = {2n + 1)7tT and Pg = 2mTT, respectively, where n is an integer. Similarly, a spatial 
Wilson loop Wx^z in real space, which involves, say the x and z direction, transforms into 
a time like Wilson loop Wj^ ^ in funny space. Hence, it can be associated with a heavy 
quark potential in funny space. The space-like Bethe-Salpeter amplitudes also obtain 
a 'physical' meaning in funny space. They can be understood as wave functions of an 
quark- ant iquark pair. 



First LGC with pure glue [|I0| and more recent calculation involving dynamical fermions 
|r3| have indicated that the space like string tension remains constant between zero tem- 
perature (where it is identical with the time like string tension) and temperatures slightly 



above the phase transition. In very recent work by Bali et al. |]T2| have studied the 



temperature dependence of the space like string tension in SU(2) gauge theory with high 
accuracy. For small temperatures, T < 2Tc, they confirm that the string tension is essen- 
tially independent of temperature. For temperatures T > 2Tc they find that the string 
tension scales like cr{T) = {g{T)T^yao where the T-dependence of the coupling g{T) 
is given by the 2-loop /5-function with a scale parameter A = 0.076Tc. This form of the 
temperature dependence is what one would expect from dimensional reduction arguments 



TM and the resulting value for the string tension agrees within 10% with that of three 



dimensional SU(2) gauge theory [|T^ . 



With the knowledge of the space like string tension as a function of temperature, the 
resulting Bethe-Salpeter amplitudes can be calculated in the model proposed in ref. p. 
This will provide a prediction for the temperature dependence of these correlations which 
can be tested on the lattice in order to check the validity of this model. In this article 
we will present the wave functions for several temperatures. We will also discuss several 
possibilities to better probe the detailed structure of the quark-antiquark interaction and 
we will investigate the implications for observables such as the screening masses etc. 

This paper is organized as follows: In the following section we briefly review the 
method we are going to use to calculate the Bethe-Salpeter amplitudes and present the 
resulting amplitudes for different temperatures. Then we investigate other possibilities 
to test the interaction more precisely. In section three we discuss the implications of the 
quark-antiquark interactions on the screening masses. In the final section we will compare 
our results with already existing LGC on the temperature dependence of the wavefunc- 



tions [|16[ and discuss possible discrepancies. 



2 Wavefunctions 

As mentioned in the introduction and explained in great detail in ref. at high tem- 
perature the effect of the space like string tension on the Bethe-Salpeter amplitude is 
best studied in 'funny space'. There the space like string tension acts just like a regu- 
lar potential and the boundary conditions, originally in the temporal direction, lead to 
standing waves in the Zf direction. At high temperatures only the lowest momentum 
(Matsubara) modes contribute, i.e. Pz = ±ttT for the quarks/antiquarks, and the mo- 
menta of the standing waves \pz\ simply act like an effective mass, as far as the motion 
in the transverse direction {xf,yf) is concerned. Under these assumptions the space like 
Bethe-Salpeter amplitude can be identified with the wavefunction of a quark-antiquark 
pair with masses m,q = ttT in two dimensions (because the motion in the Zf direction leads 
to the effective mass) interacting via the space like string potential. At high temperature, 
finally, the effective quark mass becomes large so that a nonrelativistic approximation can 
be made. 

Unless otherwise noted, from now on all variables and equations are understood to 
be defined in 'funny space'. The subscript 'f, indicating funny space variables, will, 
therefore, be dropped. 

The Schrodinger equation to be solved is M: 



1 



r dr \ dr I r^ 



ij + V{r)ij = Eip (1) 



where / is the angular momentum reflecting the azimuthal symmetry and V(r) is the 
potential which can be extracted from the space like Wilson loop. Before we will show 
the actual wavefunctions let us first discuss some general properties. For a stringlike 
potential V{r) = —g/r + err + const the large distance behavior of the wavefunction is 
governed by 

/ 1 rf2 \ 

+ ar W- = (2) 



Y ttT dr"^ 
and consequently 

^P ~ e^^^^'^' r^oG (3) 

From this asymptotic form it is evident that the wavefunctions should become narrower 
with increasing temperature even if the string tension does not increase! This is simply 
due to the increase of the effective quark mass nig = ttT in 'funny space' or, equivalently, 
to the decrease of the de Broglie wavelength of the quarks. An increase in the string 
tension would lead to an even stronger fall off of the wavefunctions. 



In fig. 1^ we show the ground state wave functions for temperatures T = .5, . . . , 8Tc 
using the temperature dependent potentials from ref. |12|, which are plotted in fig. [l|. For 



comparison in fig. ^ the wavefunctions have been calculated with a potential fixed at T = 
0.5Tc. In both cases the wave functions become narrower with increasing temperature, 
somewhat more in the first case where the string tension increases with temperature as 
well. The difference between the wavefunctions in (a) and (b) for a given temperature is 
not very large, however. As suggested by eq. (|^) a better way to illustrate the effect of 
the temperature dependent potential is to plot the wavefunctions as function of r^'^T^'^. 
This is done in fig. ^. The wavefunction in the constant potential (fig. ^) nicely follow 
the above scaling law, eq. (|^), (in case of T = 8Tc one only can observe the concave 
shape which would eventually lead to a slope parallel to the others.) In fig. ^, on the 
other hand, the effect of the string potential is clearly seen. Although in principle the 
string tension can be read off from the asymptotic slope, it requires the knowledge of the 
Bethe-Salpeter amplitudes over several orders of magnitude and thus may not be feasible 
in practical LGC As can be seen from fig. |I], over the range of the de Broglie wavelength 
d ~ 1/T, where the bulk of the wavefunction is expected to be located, the effect of 
temperature on the potential is small. Thus, in order to efficiently probe the long range 
part of the interaction, one has to find ways to move the wavefunction to larger distances. 
One such possibility is to project on states with finite orbital angular momentum (/ 7^ 
in eq. (|l|)). These states 'feel' a centrifugal barrier ~ P/r"^ and, hence, are pushed outside 
to larger distances as can be seen in fig.^, which shows the wavefunctions for 1=0,1,2,3 for 
a temperature T = 2Tc. The increased sensitivity on the potential is demonstrated in fig. 
p. There the ratio of the average radii for wavefunctions in the temperature dependent 
and temperature independent potential {< r{T — dep) > / < r(T — indep) >) with 

<r >= I d^rr\il){r)\'^ (4) 

is plotted as a function of the temperature for orbital angular momentum / = (open 
squares) and (1=3) full squares. In case of / = the average radius is almost insensitive 
to the temperature dependence of the potential. Even for the highest temperature, where 
the string tension has increased considerably, the difference in the mean radii of the wave 
functions is only ~ 10%. Thus, due to the increase of the quark effective mass, the (/ = 0) 
wavefunctions only 'see' the short distance part of the potential, where the temperature 
dependence is small (see fig. (|I])). Only the tails of the wavefunctions are sensitive to 
the linear rising part of the potential as demonstrated in figs. ^ and ^. These of course 
affect the mean radius only very little. In case of the / = 3 wave functions, on the other 
hand, the centrifugal potential /SV = /^/r^ dominates at short distances and pushes the 
wavefunctions out into the region where it becomes sensitive to the string potential. At 
temperatures T = T^ and T = 1.33 T,,, where the string tension differs very little from the 



T = 0.5Tc value, an effect in the ratio of the mean radii can aheady be observed. At high 
temperatures, the 'squeezing' of the state due to the increased string potential becomes 
very clearQ. 

In summary, the bulk of the s-wave (/ = 0) wave functions is mainly dominated by 
the effective quark mass and information about the strength of the string tension can be 
extracted only from very precise LGC The sensitivity can be increased, however, by pro- 
jecting on states with finite angular momentum. Thus sources and sinks in LGC should 
be designed in a way that they produce states of good angular momentum in the two 
dimensional transverse (x,y) plane. Similar projections have already been carried out in 
the study of the glueball spectrum. 



3 Screening masses 

In principle one would expect the temperature dependence of the space like string potential 
to reflect on the screening masses. The eigenvalues Eschr of the Schrodinger equation are 
directly related to the screening mass by 

niscr = 27rT + Eschr (5) 

and an increase of the potential results in an increase of the eigenvalue. But the eigenvalue 
is also affected by the T-dependence of the effective quark mass, because with increasing 
mass the wavefunction can sit deeper in the potential well. This dropping of the eigenvalue 
as a result of the increasing quark mass is demonstrated by the dashed line in fig. (^, 
which gives the eigenvalue as a function of temperature for the constant potential. If 
the temperature dependence of the potential is taken into account as well (full line) the 
eigenvalue does not change with temperature for T > 2 T,,. Obviously the decrease of the 
eigenvalue due to the effective mass and the increase as a result of the potential cancel 
each other with a surprising accuracy. 

But more importantly, the screening mass is dominated by the 'thermal' mass term 
~ 27rT. The energy eigenvalue Eschr is only a small correction to this termQ and the effect 
of the temperature dependence of the potential is even smaller. This is demonstrated in 
fig. ^ (lower curves), where the ratio Mscr/27rT is plotted. For all temperatures both 

^In principle the wavefunctions of the radiaUy excited states with I = also extend further out (see 
fig. 0) and, hence, should be more sensitive to the string potential. However, it seems very difficult to 
device sources and sinks which are sensitive to one particular radial excitation only. 

•^The nonrelativistic approximation, which led to the Schrodinger equation (nh, implies, of course, that 
the mass term is much larger than the nonrelativistic energy eigenvalue. 



temperature dependent and temperature independent potential lead to the same screen- 
ing mass. Consequently, any effects from the T-dependence of the string potential on the 
screening mass, although possibly present, are not accessible on lattices of the present 
size. For example for T = AT^, AEschr — 0.5 Tc and Mgcr — '^5Tc. Hence, in order to 
be sensitive to the temperature dependence of the potential a lattice which is about 200 
times as large in the spatial direction than in the time direction would be required ! 



4 Point sources 

Usually screening masses are defined as the logarithmic derivative of point to point corre- 
lators taken at large distances. On the lattice, however, the maximum separation between 
source and sink of the correlator is limited for practical reasons. In this case the signal 
may be 'contaminated' by an admixture of excited states, because point source have an 
overlap with all radial excitations of the quark antiquark-system and the limited separa- 
tion between sink and source may not be enough to screen the excited states sufficiently. 
As we will show, in our case, the presence of excited states hardly affects the results for 
the screening masses but leads to considerable corrections for the Bethe-Salpeter Ampli- 
tudes. Assuming the validity of the above potential model, the point to point correlation 
function in funny space can be written as (for details see appendix) 

C(r) = -^El^n(0)Pe-^'"^ (6) 

with M„ = 27rT + Eschr and V4 is the volume of the four-dimensional euclidian box. 
Here the contribution of states with higher Matsubara frequencies has been neglected 
(see appendix). Because of the factor |'?/'(0)p only those states with angular momentum 
/ = contribute to the sum in eq. (^. Similarly, for the Bethe-Salpeter amplitudes one 
has (see appendix) 

^{r,r^) = ^Y.r{Q)i'{ri.)e-^-^ (7) 

As already mentioned, the ground state is projected out by going to large distances r. 
Usually this distance is assumed to be reached once the screening mass 

M,,, = ^ln(C(r)) (8) 

does not vary with r anymore. We will demonstrate that this method may have practical 
problems, if one wants to extract the ground state wavefunctions as well. The argument 



is, that, due to the dominance of the ' free, thermal' contribution ~ vrT to the screening 
mass, the screening mass of the ground state and of the first excited state differ only by a 
few percent. Therefore, the screening mass appears to have assumed its asymptotic value 
even if a nonnegligible admixture of excited states is still present. While these excited 
states obviously do not alter the value of the screening mass, they lead to substantial cor- 
rections for the Bethe-Salpeter amplitudes. The reason is, that at intermediate transverse 
distances r_|_ the contributions to the 'wavef unction' of the ground- and first excited state 
(n = 1, / = 0) are comparable in magnitude but opposite in sign. Thus, the resulting 
wavefunction appears narrower than that of the ground stateQ. This is demonstrated in 
fig. 1^, where in part (a) we have plotted the relative difference of the local screening 
mass Mscr(yT) (eq. 4.3) to that of the true ground state Mscr{oo) as a function of r for 
temperatures T = Tc and T = 1.33Tc0. (For the results presented here, the four lowest 
states (n=0,l,2,3) have been taken into account.) At r = l/T^ the resulting screening 
mass already differs by less than 3 % from that of the ground stateQ. In (b) the same 
ratio is plotted for the average radii of the wavefunction as defined in eq. (^). At the 
same distance r = 1/Tc the difference is still of the order of 20 %. This difference in 
the wavefunction is further illustrated in fig. || where we have plotted the radial density 
distributions r±\ip{r±)\'^ for the true ground state wavefunction (thick lines) and for the 
wavefunction as seen at a distance r = 1/Tc. The narrowing of the wavefunction can 
clearly be seen. This narrowing actually has been observed in the work of ref. [jl6[ where 



the wavefunction as obtained from a LGC, using point sources, is compared with the 
ground state wave function of the potential model for temperatures T = 0.66, 1.0, 1.33Tc. 
The lattice used has a spatial extent of exactly 2/Tc, resulting in a maximum distance of 
Tf = Tc in funny space, because of periodicity. Thus, the conditions are identical to those 
considered here. For the highest temperature, where the approximations of our potential 
model are best justified, agreement could only be found if one assumed a quark mass 
which is twice as big as given in the modelQ. Although a fairly large current quark mass 
rriq = 200 MeV has been used, this cannot account for doubling the effective quark mass 

^The narrowing of the wavefunction due to the admixture of excited state is simply a remnant of 
the initial (5-function, which corresponds to a sum over all excited states, due to completeness. With 
increasing t this ^-function disperses and eventually assumes the shape of the ground state. 

^We have selected these comparatively low temperatures, because later on we want to discuss the 
results of U& , where the same temperatures have been chosen. The choice of temperature does not affect 
our argument and we find the same results at higher temperatures. 

^Usually in LGC periodic boundary conditions are imposed in the real-space spatial directions, which 
implies that only half the lattice can be used to relax the wavefunction in funny space. Therefore, 
distances r much larger the l/T^ are rarely realized. 

^In ref. ||l6| the calculation in the potential model accidentally has been made using the full quark 
mass in the Schrodinger equation instead of the reduced mass |17|] . Thus, the agreement the authors find 
for the highest temperature in effect means that one has to double the quark mass. 



of ttT at all. The fact that a large effective quark mass is needed in order to reproduce 
the lattice result implies, that the lattice wavefunction is narrower than the one given by 
the potential model. This difference, of course, is exactly what we would predict if an 
admixture of the excited states is still present. For comparison in fig. ^ we have also 
plotted the ground state wave function obtained after doubling the quark mass (dashed- 
dotted line). We find that it agrees fairly well with the effective wavefunction extracted 
at a distance r = 1/Tc (thin dashed line). 

At this point we should note, that the narrowing of the apparent wavefunction due to 
the superposition of excited states may also be relevant for the measurement of Bethe- 
Salpeter amplitudes at zero temperature. Monitoring the screening mass in order to 
assure sufficient screening of the excited states is only successful if the energy eigenvalue 
of the ground state is smaller than the excitation energy of the first radial excitation 
(see appendix for a more formal argument). This implies that the measurement of the 
Bethe-Salpeter amplitudes for the pion using this method should be fine while in case of 
the p the monitoring of the wavefunction itself is advised. 

The authors of ref. |jl6| also find that the wavefunctions obtained on the lattice for the 



three different temperatures (T = 0.66, 1, 1.33Tc) are almost identical. In contrast our 
model would predict a narrowing of the wavefunctions due to the increase of the effective 
quark mass ~ ttT (see fig. |]). (The changes in the potential are minimal over the 
temperature interval considered here.) Also for the more realistic wavefunctions, where 
the presence of excited states due to point sources is taken into account, the narrowing as a 
function of the temperature persists. The ratio of the average radii between wavefunctions 
of different temperatures is about the same for the true ground state wavefunction and 
that taken at a distance r = l/T^ (see also fig. |^). If we assume that the qualitative 
behavior of the space like string potential is the same in SU(3) and in SU(2) gauge 
theory, the independence of the wavefunction on the temperature as observed in ref. 
||T6| cannot be understood within our model. Of course the temperatures considered 
T < 1.33Tc are somewhat low for the dimensional reductions arguments of our model to 
be good. Furthermore, at the lower two temperature points, T = .66Tc and T = T^, the 
constituent quark mass may still be of considerable magnitude, which also would affect the 
wavefunction (the authors give a screening mass for the p at T = .66Tc of about ~ 5.7Tc 
which should correspond to roughly twice the constituent quark mass). The results of ref. 
T6| , therefore, may be at or below the lower limit of applicability of our model. However, 



it would be very interesting to carry out LGC calculations at larger temperatures T > 2Tc 
in order to see if the wavefunctions eventually follow the prediction of our model or if the 
independence of the temperature, as found in ref. [|16| persist even at temperatures where 
one would expect dimensional reduction to be valid. 



5 Conclusions 

In this article we have studied the effect of the temperature dependence of the space-hke 
Wilson loop on the Bethe-Salpeter amplitudes in the space-like direction. The tempera- 
ture dependence of the wavefunctions is controlled not only by the temperature depen- 
dence of the potential but also by the increase of the effective quark mass with temperature 
in funny space {nig = irT). Because of that 'trivial' effect, only high accuracy LGC will be 
able to distinguish between ground state wavefunctions from temperature dependent and 
temperature independent string potential at high temperatures. We, therefore, propose 
to project instead on states with finite angular momentum in the transverse directions. 
As a result of the centrifugal potential, these states are pushed towards larger transverse 
distances and are thus much more sensitive to the string potential. 

We have further demonstrated, that screening masses are extremely insensitive to the 
temperature dependence of the potential. They are essentially dominated by the 'free 
thermal' value ~ 27rT and deviations due to the string potential are very small. We 
also have shown that the presence of excited states in the correlation function due to 
point sources minimally changes the result for the screening mass. The Bethe-Salpeter 
amplitudes, on the other hand, come out much more narrow if excited states are still 
present. As a consequence, monitoring the screening mass in order to decide if a state at 
a given distance has settled into the ground state is not sufficient and may lead to wrong 
conclusions about the actual size of the Bethe-Salpeter amplitudes. 

Finally, we have presented unique predictions for the temperature dependence of the 
space like Bethe-Salpeter amphtudes based on the potential extracted from high accuracy 
SU(2) LGC. In our model ^ these amplitudes are understood as Schrodinger wavefunc- 
tions of quarks with mass vrT confined by the space like string potential. A measurement 
of these amplitudes on the lattice would, therefore, provides an important test for the 
validity of our understanding of these correlations. 
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6 Appendix 

In this appendix we will derive how point to point correlation functions and Bethe-Salpeter 
amplitudes are expressed in terms of the eigenstates and eigenvalues of the Schrodinger 
equation (|l|). We will also give an estimate on how the presence of excited states affects 
the results for screening masses and Bethe-Salpeter amplitudes. In the entire appendix 
we will work in funny space i.e. we have periodic/antiperiodic boundary conditions in 
z-direction and none in the r-direction. 



6.1 Point to point correlation functions 

The screening masses are determined from the exponential fall off in the r direction of 
the point to point correlation functions. The correlation function is given by: 

C(r,f)=<0|j(r,r-')j(0)|0> (9) 

where r > and r^ = z > is assumed. Our main assumption is that this correlation 
function can be saturated by the eigenstates of the Schrodinger eq. ([l|). Therefore, we 
can insert a complete set of states 

where p refers to the momentum of the center of mass motion and n to the quantum 
numbers of the internal motion, i.e. to that of the Schrodinger equation. The correlation 
function can then be written as 

Cir,r) = Yl TT^T < 0\j{T,f)\n,p X n,p1j(0)|0 > (11) 

Since \n,p > are eigenstates of the momentum operator - with discrete eigenvalues for 

pz- we can write 

C{T,f) = E l< 0\j{0)\n,p>fexp{-En{p)T + tpf) (12) 

n,p 

By integrating over the spatial directions we project on states with zero transverse mo- 
mentum and zero Matsubara frequency Pz of the quark antiquark pair. 

C{r) = Y. |<0|j(0)|n,p = 0>|'exp(-E„r) (13) 
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The sum still includes all pairs of quark Matsubara frequencies which add to zero, i.e. 
((2n + 1)itT; — (2n+ l)7rT). We will, however, neglect the contribution of the states with 
\n\ > 1, because they have a screening mass of Mg^r > 47rT and, thus, are screened very 
efficiently. 

Finally the matrix element < 0\j{0)\n,p> relates to the Schrodinger wavefunction by 



< 0|j(0)|n,p >= ^2E4p)^^,p{0) (14) 

with 

vl>(f) = -L$„_^(f)^„(o) (15) 

Here $(r) denotes the wavefunction of the center of mass motion and ipn{0) is the wave- 
function of eigenstates of the Schrodinger equation describing the relative motion in the 
transverse (x, y) direction. (3 = 1/T gives the box size in the z-direction in funny space. 
The normalization factor 1/v^ is the remnant of the standing waves ~ exp(-tp^z) ^^^ ^^^^ 

motion in the z-direction in the relative coordinate. Therefore, with |$„(0)p = I/V3 the 
correlation function is given by (V4 = V3/3) 

C{r) = ^Y: \MO)fexp{-E^T) (16) 

with En = 2ttT + Eschr, where Eschr is the eigenvalue of the Schrodinger equation ([|). As 
explained above, only the lowest Matsubara frequency of the individual quarks are taken 
into account. 

6.2 Bethe-Salpeter amplitudes 

The Bethe-Salpeter amplitude is defined as 

5(r,f,rx)=<0|(^(r,r-^)r^(r,r + ^)) j(0)|0 > (17) 

where F denotes the appropriate spin-isospin matrix, which we will, for simplicity, ignore 
from now on. Again we can write the product 



<rz,p1^(r,r-^)F^(r,r + ^)|0>= ^ ^^^^ $„,,(rXr^) (18) 

(Note, that r±_ is in the transverse (x, y) direction only.) By inserting a complete set of 
states, we can factor out the cm. motion and obtain 

B{T,f,r^) =Y.< n,p\m\0 > '-^^^^^^=^^Mr±) (19) 

n,p V^-t^n 
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and after integrating over the spatial coordinates the states |n,p = > with vanishing 



cm. momentum are projected out. Using eq. (|T^ we finally have 



B(T,r^) 



Va 



E<(0)V'nK)exp(-E„r) 



(20) 



6.3 Effect of excited states on screening mass and Bethe-Salpeter 
amplitudes 

Finally let us give an estimate on how the presence of excited states affects the screening 
mass and the Bethe-Salpeter amplitude. Let us, for simplicity, consider only 2-states. 
Then the correlation function is given by (up to a constant) 



C(r) = exp{—EQT) + a^exp(— i^ir) 



with 



a 



MO) 



MO) 



(21) 



(22) 



The screening mass at given distance r follows from the logarithmic derivative of the 
correlation function 



with 



M. 



A = a 



ar 



2 ^1 - Eq 
Eq 



exp(-(Ei-Eo)r) 



(23) 



(24) 



The screening mass is, thus, unaffected by the presence of the excited state if |A| ^ 1. 
This is achieved by going to a sufficiently large distance r. 

For the Bethe-Salpeter amplitude or wavefunction, on the other hand, one has 



?/'(r_L) 



a 



Mr±) exp(-Eor) + a^i(r_L) exp(-Eir) 
MO) 



± 



V^i(o) 



(25) 



which can be related to the accuracy of the screening mass A 

A En 



V'(^±) = Mr±) exp{-EoT) ( 1 ± — - — ^ ) 

V a"^ El — EfjJ 



(26) 
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Typically a^ < 1 (see e.g. fig. ^ and hence the accuracy of the Bethe-Salpeter amplitude 
As.s, is 

Ab.s. > As.M. ^ ^° (27) 

Therefore, the error in the wavefunction may be considerably larger than that in the 
screening mass if the energy eigenvalue of the ground state is much larger than the exci- 
tation energy of the first excited state, i.e. Eq ^ Ei — Eq. 
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Figure captions 



Figure 1: Space like string potential from ref.|12 



Figure 2: Space like Bethe-Salpeter amplitude for ground state for temperature dependent 
potential (a) and for potential fixed at T = 0.5 Tc (b). 



Figure 3: Space like Bethe-Salpeter amplitude for ground state for temperature dependent 
potential (a) and for potential fixed at T = 0.5 Tc (b). Here x-axis has been rescaled to 
demonstrate the asymptotic behavior. 
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Figure 4: Space like Bethe-Salpeter amplitude for ground state (n=0,l=0), the first radial 
excitation (n=l,l=0), and for finite angular momentum (n=0,l=l,2,3). The temperature 
is T = 2 T,. 



Figure 5: Ratio of average radii between temperature dependent and temperature inde- 
pendent potential for states with 1=0 (full line) and 1=3 (dashed line) 



Figure 6: Eigenvalues of the ground state of the Schrodinger equation and associated 
screening masses for temperature dependent (solid line) and temperature independent 
(dashed line). The dotted line gives the screening mass of the first radial excitation. 



Figure 7: Apparent screening mass (a) and average radii of the apparent wavefunctions 
as a function of separation of source and sink r 



Figure 8: Probability density distribution for the true ground state (thin lines) and for 
the effective wavefunction as obtained at a distance of r = 1/Tc (thick lines). The 
dashed-dotted line corresponds to the wavefunction obtained with twice the quark mass 
for T = 1.33 T, 
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